Magnetostriction and exchange effects in trapped dipolar Bose and Fermi gases 
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We examine the magnetostrictive position and momentum space distortions that occur in harmonically con- 
fined dipolar Bose and Fermi gases. Direct interactions give rise to position space magnetostriction and ex- 
change interactions give rise to momentum space magnetostriction. While the position space magnetostriction 
is similar in Bose and Fermi systems, the momentum space magnetostriction is markedly different: the Bose gas 
momentum distribution distorts in the opposite sense to that of the Fermi gas. By relating exchange effects to 
short range correlations between the particles we discuss the energetic origin of this difference. Our main calcu- 
lations are based on Hartree-Fock theory, but we also provide analytic approximations for the magnetostriction 
effects at zero and finite temperature. Our predictions should be verifiable in current experiments with ultra-cold 
polar molecules. 
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I. INTRODUCTION 

Recently a range of ultra-cold gases with appreciable mag- 
netic [1^1] and electric [5, 6] dipoles have been realized in 
experiments. The key feature of these gases is that the con- 
stituent particles interact via a dipole-dipole interaction (DDI) 
that is long-ranged and anisotropic. When these systems are 
polarized by an external field the DDI causes their shape 
to distort. The first experimental evidence for such a mag- 
netostrictive effect 1 was in the expansion of dipolar Bose- 
Einstein condensates (BECs) [7] (also see [2, 8-10]). 

Here we examine magnetostrictive effects in normal Bose 
and Fermi gases with polarized dipoles. In our work we dis- 
tinguish between two cases: position space magnetostriction 
and momentum space magnetostriction, whereby the posi- 
tion and momentum space distributions, respectively, distort 
relative to the direction of the polarized dipoles. We show 
that position space magnetostriction is qualitatively similar for 
Bose and Fermi gases: the gas tends to elongate along the 
dipole polarization direction to reduce the direct interaction 
energy. Momentum space magnetostriction arises solely from 
exchange interactions and distinguishes the quantum statistics 
of the particles: the Fermi gas momentum distribution elon- 
gates along the polarization direction, while the Bose gas in- 
stead contracts along this direction. 

Previous studies have considered the distortion of ultra-cold 
gases with DDIs. Work on Bose gases has focused on posi- 
tion space magnetostriction for zero temperature BECs [11- 
13]. For normal Fermi gases Hartree calculations (excluding 
exchange) were used to study position space magnetostric- 
tion in [14]. The first theories including exchange interactions 
for the trapped Fermi gas were based on simple variational 
approximations [15-18] (also see [19]), and revealed distor- 
tion of the momentum distribution. In this paper we use full 
Hartree-Fock (HF) calculations to provide a unified quantita- 
tive treatment of trapped Bose and Fermi gases. While HF 
calculations have been applied to trapped Fermi gases [20- 



22], this is the first application of this theory the trapped dipo- 
lar Bose gas (c.f. Hartree analysis reported in [23], also see 
[24]). Our main concern is in characterizing magnetostrictive 
effects, their relationship to direct and exchange interactions, 
and to pair correlations. We support our numerical results by 
analytical approximations. 



II. THEORY 

We consider a gas of polarized dipolar particles that interact 
by a long-range dipole-dipole interaction of the form 



C/dd(x) 
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where the constant Cdd is given by ^oMm f° r magnetic dipoles 
of strength \i m and d 2 /e for electric dipoles of strength d, and 
9 is the angle between x and the polarization axis, which we 
take to be the z direction. The atoms are confined within a 
cylindrically symmetric harmonic trap 
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with aspect ratio X = uj z /uj p . 

Following [20-22, 25] we introduce the dimensionless in- 
teraction parameter D t = CmN 1 ^ 6 / ' {Airhuia\ ) where w = 
(ujpLUz) 1 / 3 and aho = \fhpmuj. When considering ther- 
mal effects we characterize the temperature either in terms of 
the ideal condensation temperature, Tj? = y/N/((3)hw/kB, 
with ((s) the zeta function, or the Fermi temperature, Tp = 
^6Nhu}/kB- We note that these two degeneracy tempera- 
tures relate as T c ° w 0.52T F . 



A. Semiclassical theory 

In the semi-classical approximation [14, 21, 22] the system 
is described by the Wigner function 



1 We adopt terminology appropriate to magnetic dipole gases. We note that 
for polar molecule gases the effects we predict are electrostrictive. 
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where 77 = 1 for bosons and rj = — I for fermions, and ll is the 
chemical potential. This approximation furnishes a good de- 
scription when the temperature is large compared to the trap 
level spacing (fcgT » huj). For fermions this approxima- 
tion can be applied at T = as long as the Fermi energy is 
sufficiently large (i.e. fi » huj) [20]. For the case of bosons 
at T < T c , where T c is the condensation temperature [26], 
a condensate emerges and the semi-classical theory requires 
augmentation by a beyond-semi-classical description of the 
condensate. Here we restrict our attention to normal Bose and 
Fermi gases. 

The position space and momentum space densities are 
given by 



n(x) 
n(k) 



f—w( 
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respectively and the total number is given by N — J dxn(x). 



the direct and exchange interaction terms, respectively. In 
Eq. (10), t/dd(k) = Cdd(cos 2 #k — 1/3) is the Fourier trans- 
form of the dipole-dipole interaction, where # k is the angle 
between k and k z . The HF theory requires solving Eqs. (3), 
(8)-(10) in an iterative procedure until the equations are self- 
consistent. We are interested in solutions for a specified total 
number, N, which requires an additional iterative procedure 
to adjust p. 

Our system is cylindrically symmetric, since the symmetry 
axis of the harmonic trap coincides with the polarization di- 
rection. We can exploit this symmetry to simplify the Wigner 
function to a function of four variables {p, z, k p , k z }, i.e., 
cylindrical coordinates in x and k space. Our numerical algo- 
rithm makes use of Fourier transform techniques and quadra- 
ture based on discrete cosine and Bessel functions, and is de- 
tailed in [22] (also see [21]). The numerical implementation 
of HF theory requires extensive computational resources be- 
cause of the necessity to construct the full Wigner function on 
grids with suitable resolution to accurately evaluate the DDI. 



B. Characterization of magnetostriction 



D. Hartree theory 



Our primary concern is distortions in the trapped clouds 
arising from the DDI. It is convenient to characterize this in 
terms of the rms-widths 

-1 1/2 
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with v = {x, y, z, k x , k y , k z }. 

In terms of these quantities we define the momentum and 
position distortion as 

1 <Jx 
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respectively, i.e. as the ratio of the widths in the directions per- 
pendicular (x or k x ) and parallel (z or k z ) to the dipolar polar- 
ization direction. The trap anisotropy appears in the definition 
of j3 so that in the absence of interactions, where the posi- 
tion density has the anisotropy imposed by the trap, we have 
j3 = 1. In contrast the non-interacting momentum distribution 
is spherically symmetric. We refer to systems as exhibiting 
position space magnetostriction when j3 ^ 1 and momentum 
space magnetostriction when a^l. 



We note that the direct interaction potential (9) can be eval- 
uated in terms of the position space density, i.e. $_d(x) = 
/ dx' t/dd(x— x')n(x'), thus if the exchange term is neglected 
the k dependence of the dispersion (8) is quadratic and Eq. (4) 
can be evaluated [22, 23] 

n( X ) = ^-Q /2 ( e [M-V «(x)]/* B T) (n) 

where C2( z ) = Yl'jLi rf -1 ^ /j a is the polylogarithm func- 
tion, A d B = hj \J2-Kmk B T , and 



V eff (x) = C/ tr (x) + $ D (x), 



(12) 



is the effective potential. By neglecting exchange we have ar- 
rived at a Hartree theory for the system, which involves solv- 
ing Eqs. (11) and (12) self-consistently. This theory is compu- 
tationally much simpler because it does not require the evalu- 
ation of the full Wigner function. The Hartree theory provides 
a baseline for investigating the exchange effects included in 
the full HF calculations. 



C. Hartree-Fock theory 



E. Analytic approximations 



The HF dispersion relation is [21, 22, 27] 
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e(x, k) = — + f/ tr (x) + $ D (x) + 77$ B (x, k), (8) 



with 
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To develop an analytic approximations for the magne- 
tostrictive behavior we introduce the following ansatz for the 
HF dispersion relation: 



e^(x,k) = ^- ( Kp k 2 p + K z k 2 z ) + — ± ( Xp p 2 + X,A 2 ~~ 1 
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where K p , k z ,x p an d Xz are variational parameters to be de- 
termined. 
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1. Zero temperature fermions 

We first consider an analytic treatment of T = fermions. 
For this system all single particle states up to the Fermi energy 
are occupied. We use a Thomas-Fermi ansatz for the Wigner 
function (after [15]) 

Wb(x,k) = e[/i-c A (x,k)], (14) 

where 0[ ] is the Heaviside step function. In the appendix, we 
minimize the energy to determine the variational parameters, 
from which we obtain a and /?. We find, to first order, that 

awl-c a D tl (15) 
/3« l + c p \- 2 J'(\- 2 - 1)A, (16) 

o3S/3 r>35/3 

where c a = 3l7 % 175 ^ 2 ~ 0.167, c p = 3 if /6357T 2 ~ 1-26, and 
the function J(u) is defined in Eq. (A3). 

2. Finite temperature 

To account for thermal effects in Bose and Fermi gases we 
use a Boltzmann ansatz for the Wigner function (after [19, 
28]) 

W T {X, k) = e [M-^(x,k)]/fc B T_ 

This approximation fails to account for the quantum statistics 
(e.g. phase space densities greater than unity are not prohib- 
ited for Fermions), and thus is only applicable at temperatures 
well above the relevant degeneracy temperatures. 2 In the ap- 
pendix, we minimize the free energy 3 derived from (17) to 
determine the variational parameters and hence the distortion 
parameters. To first order we find 

,. 1+ ^ w - 1) (^) 5 % t . m 
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FIG. 1. (Color online) Position space distortion for fermions using 
HF (solid curves) and Hartree (diamonds), and for bosons using HF 
(dashed curves) and Hartree (circles), also the predictions from (16) 
(crosses) and (19) (thin curves) with A = 0. 1 [blue (dark gray)], 1 
[green (light gray)], and 10 [red (gray)], (a) Bosons and fermions at 
T = 1.5T C ° w 0.78T£. (b) Fermions with D t = 1. (inset) Fermions 
at T = 0. 



III. RESULTS 

A. Position space magnetostriction 

In Fig. 1(a) we compare the HF position space distor- 
tion against the analytic expression (19), and Hartree results. 
For all results, f3 decreases with increasing dipole interaction 



2 We are unable to develop simple closed-form analytic results using de- 
generate forms of the Wigner function [c.f. Eq. (3)]. Nevertheless, such 
a degenerate form could be treated numerically and would lead to better 
agreement with the full HF results. 

3 We note that Bose and Fermi statistics are included in the treatment of 
interactions by adding the appropriate exchange interaction term. 



strength, i.e. the gas distorts so as to extend along the direc- 
tion that the dipoles are polarized. This distortion increases 
the number of dipoles in the energetically favorable head-to- 
tail configuration. The numerical results terminate at finite 
DDI due to the system nearing mechanical instability (e.g. see 
[23]). We observe that at any given value of D t the position 
distortion is greatest in near-spherical geometries. The ana- 
lytic results are independent of whether particles are bosons 
or fermions, and predict a linear change in j3 with the interac- 
tion strength. The Hartree results are in good agreement with 
the full HF solutions for the distortion. This demonstrates that 
the position space distortion arises from the direct interaction. 

The T~ 5 / 2 scaling in the finite temperature analytic pre- 
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diction (19) suggests that large magnetostriction effects will 
be obtained at low temperature. For bosons the critical tem- 
perature for condensation T c forms a lower bound for the ap- 
plicability of our theory. For fermions our theory could be 
applied as T — > 0, however the non-degenerate form of the 
distribution used to derive the analytic expression (19) will 
be unreliable below the Fermi temperature. In Fig. 1(b) we 
compare the finite temperature analytic results to the HF cal- 
culations for a range of temperatures below the Fermi temper- 
ature. These results indicate that the analytic approximation 
(19) significantly overestimates the distortion for T < 0.5T F , 
and diverges as T — > 0. We emphasize that this breakdown is 
not entirely due to the first order expansion used to derive the 
analytic results, but also because the distribution (17) fails to 
work in the degenerate regime. The T = analytic expres- 
sion (16) is also shown in Fig. 1(b) and provides a qualitative 
estimate of the HF calculations at T — >• 0. However, for the 
interaction strength used in Fig. 1(b) (D t — 1) corrections to 
our linearized result are apparent. In the inset to Fig. 1(b) we 
compare the HF results to the analytic approximation at T = 
showing that good agreement is obtained for smaller values of 
D t . 



B. Momentum space magnetostriction 

In Fig. 2(a) we compare the full HF results for the momen- 
tum space distortion against the analytic expression (18). For 
the fermionic system a decreases with increasing dipole in- 
teraction strength, i.e. the gas distorts in momentum space to 
increase its extent in the k z direction. In contrast, for the 
bosonic system, a increases with increasing dipole interac- 
tion strength, thus the system contracts along the k z direction 
relative to its radial extent. The Hartree theory (not shown) 
predicts no momentum distortion (i.e. a — 1), thus momen- 
tum distortion arises solely from the exchange interactions, 
which are neglected in the Hartree approximation. As in the 
previous section, the numerical results terminate at finite DDI 
due to the system nearing mechanical instability. 

In Fig. 2(b) we consider the Fermi system in the low tem- 
perature regime (T < T F ). Similar to the position space 
distortion [see Fig. 1(b)] we observe that the analytic ap- 
proximation (18) significantly overestimates the distortion for 
T < 0.5T£, and diverges as T -> 0. The T = analytic 
expression (15) provides a qualitative estimate of the HF cal- 
culations at T — > 0. In the inset we compare the HF results 
to the analytic approximation at T = showing that good 
agreement is obtained for smaller values of D t . 



C. Exchange interactions and second order correlations 

To understand the role of exchange effects it is useful to 
consider the second order correlation function 

G (2) (x,x') = ($t(x)$ t (x')$(x')$(x)\ , (20) 

which describes the probability of finding particles at posi- 
tions x and x', with $>(x) being the field operator. It is conve- 




FIG. 2. (Color online) Momentum distortion for fermions (solid 
curves) and bosons (dashed curves) and the predictions from (15) 
(black cross) and (18) (thin black curves) with A = 0.1 [blue (dark 
gray)], 1 [green (light gray)], and 10 [red (gray)], (a) Bosons and 
fermions at T = 1.5T C ° « 0.78T£, (b) Fermions with D t = 1, 
(inset) Fermions at T = 0. 



nient to transform this to center of mass R = ~(x + x') and 
relative r = x — x' coordinates, i.e. (R, r), and introduce 
the pair correlation function 



G(R, r) = G (2) (R, r) - n(R + |r)n(R - 



(21) 



Within HF theory, higher order correlation functions decom- 
pose into products of first order correlation functions (e.g. see 
[22, 29]) giving 



G(R, r) = 7] 



W(R,k)e 



ik-r 



(22) 



The T) prefactor of this result shows the important effect of 
quantum statistics on correlations: For fermions G(R, r) < 
indicates their tendency to anti-bunch, while for bosons 
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C(R, r) > indicates their tendency to bunch [30]. Around 
any given point R the pair correlations decay as |r| increases, 
with this decay length defining the correlation length. Equa- 
tion (22) shows that the r dependence of the correlation func- 
tion is related to the Fourier transform of the momentum dis- 
tribution. Thus anisotropy in the momentum dependence of 
W(R, k) will be accompanied by an inverse anisotropy in the 
r dependence of the pair correlation function. 

We present a comparison of C(R, r) for a Bose and Fermi 
gas in Fig. 3 that demonstrates the features discussed above. 
The inset to that figure shows the density distribution, and re- 
veals that the correlation length scale (~ AdB) is much smaller 
than the spatial extent of the systems. 

In terms of the density and pair correlation functions we 
can evaluate the interaction energy densities 

£ D (R) = ^[ drU dd {r)n(R+ §r)n(R- |r), (23) 
£ E (R) = ^J drC/ dd (r)C(R,r), (24) 

for the direct and exchange terms, respectively. Note that the 
total direct and exchange energies, obtained by integrating the 
respective energy densities, are given by 

E D = J dR£ D {R) = \J ^<f> D ( x )iy(x,k), (25) 



E E = / dR£ E (R) 



V J ^<Mx,k)iy(x,k). (26) 



Equation (23) shows that the direct energy arises from the den- 
sity distribution of the system, while the exchange term (24) 
arises from correlations, i.e. with the integrand in Eq. (24) 
only being non-zero over the small range of r values for which 
correlations extend [e.g. see Fig. 3]. 

The anisotropy of the correlations ensures that the exchange 
energy is minimized. To see this we note that the angular part 
of /7dd(r) is proportional to the Y^iO, <fi) spherical harmonic, 
so that the integral in Eq. (24) is sensitive to a quadrupolar mo- 
ment (i.e. distortion) of C(R, r) in the r variable. For bosons, 
where the correlation is positive, the correlation peak is pro- 
late: short range correlations increase the number of pairs 
in the energetically favorable head-to-tail configuration. In 
fermions, where the correlation is negative, the correlation dip 
is oblate: short range correlations reduce the number of pairs 
in the energetically unfavorable side-by-side configuration. 

Finally, while our focus here has been on second order cor- 
relations, because of their relationship to energy, the exchange 
effects are also apparent in the first order coherence function 
G^^x, x') = (x')'0(x)) which, in center of mass and 
relative coordinates, relates to the distribution as 

G«(R,r) = |^^(R,k)e lk - r . (27) 

This function also exhibits the anisotropy from exchange in- 
teractions with a differing sense of distortion for Bose and 
Fermi gases. 
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FIG. 3. (Color online) Second order correlation function for A = 10, 
D t = 4, N = 10 5 , T = 1.1T° « 0.57T£ for Fermi (solid curve) 
and Bose (dashed curve) for r = pp (blue, dark grey) r = zz (red, 
light grey). Inset shows the density profiles for comparison. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper we have examined magnetostrictive effects in 
normal Bose and Fermi gases with dipole-dipole interactions. 
We have performed Hartree-Fock calculations for these sys- 
tems and compared the results against Hartree calculations 
and analytic approximations that we have developed. Mag- 
netostrictive effects are only apparent for moderately large 
dipole strengths (D t of the order of unity), as being realized in 
current experiments with polar molecules. E.g. N = 10 5 KRb 
molecules (for which both Bose [5] and Fermi [6] molecules 
have been realized) with d = 0.5 D in a w = 2-7T x 100 s -1 
trap corresponds to D t = 4.7 (higher interaction strengths 
are possible with more polar molecules, e.g. NaK which has 
d « 2.72 D [31]). We note that Hartree-Fock theory should 
provide an accurate description of the system in the weak- 
coupling regime, e.g. in the regime of dipolar atom (mag- 
netic dipole) experiments. For polar molecule experiments, 
because of the larger dipole strengths, corrections to HF the- 
ory may become important, however it is expected that the 
HF results should be qualitatively accurate (e.g. see discus- 
sion in Ref. [32], also see [33]). Here we have focused on 
the case of pure DDIs to allow a more direct comparison of 
Bose and Fermi systems. However, we make a few general 
comments about the effect of contact interactions on the Bose 
system based on calculations we have performed. First, re- 
pulsive contact interactions allow the normal Bose gas to be 
stable out to higher dipole strengths (e.g. see [23]), where the 
DDI induced distortions are larger. Second, for a given case 
(i.e. number of particles, DDI strength and temperature) the 
addition of repulsive contact interactions tends to lower the 
density and this in turn causes the DDI induced distortions to 
slightly decrease. 

It should be possible to experimentally verify the effects we 
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have studied in this paper. Spatial distributions have been pre- 
cisely measured using absorption imaging in cold atom ex- 
periments [34, 35], and the extension of these techniques to 
the dipolar gas should allow the position space magnetostric- 
tive behavior to be characterized. The direct measurement of 
the momentum space effect is possibly more challenging. Re- 
leasing the system from the trap requires a model of the ex- 
pansion dynamics (e.g. see [17, 25, 28]) and does not simply 
reveal the momentum behavior. Bragg spectroscopy could be 
used to probe the momentum distribution in situ [36, 37]. As 
discussed in this paper the momentum space magnetostrictive 
effects in dipolar gases are revealed by the development of 
anisotropic correlation functions. These could be revealed by 
a measurement of the second order (density-density) correla- 
tion function, or by measurements of the first order coherence 
function [c.f. Eq. (27)]. Such measurements have been per- 
formed with high accuracy on alkali Bose gases [34, 35, 38]. 
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Appendix A: Analytic approximations 



1. Zero temperature fermions 



Using (14), the energy of a dipolar gas described by Wq is 
(after [15]), ignoring factors of XpXz and n 2 p n z which do not 
contribute to the distortions to first order 
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where S — k z /k p — 1 and' 
1 



4 



J(u) = 



d6 sinf 



1 + u sin 

-l 



-sinh yju 



- 1 - J(S) 



For very small values of u, the series expansion 



8v^ 
105 



I6u 6 
~315~ 



+ o(« 4 ) 



(Al) 

(A2) 
(A3) 

(A4) 



and its derivative are useful. We plot J(u) in Fig. 4(a) and 
note that it is monotonically decreasing and that it is zero in 
the spherical case, i.e. J(0) = 0. 




FIG. 4. (a) The function J(u) (solid curve) showing J(— 1) = | 
(cross), J(0) = (plus) J(u) — > — ~ as u — > oc (dashed line), (b) 
The term appearing in (16) and (19) as a function of the aspect ratio 
showing J'(0) = at A = 1 (plus). 



To find a = \/l + 5 and j3 = \JxzlXp to m " st order, we 
minimize the energy (Al), assuming 5 and Xz/Xp ~ 1 sire 
small parameters, we find Eqs. (15) and (16). We plot the 
direct interaction derivative term, A~ 2 J'(A~ 2 — 1), appearing 
in Eq. (16) [also in Eq. (19)], in Fig. 4(b) and note that the rate 
of change of direct interaction energy is largest for spherical to 
slightly oblate geometries, so that these are most able to offset 
the increase in trap energy and are likely to distort most. 

The derivative of J{u) is required to evaluate Eqs. (16) and 
(19) and is given explicitly by 



4 We note that (sinh 1 •/u) / sfu is real for u > 
+ «) 1 / 3 ]/6of [28] 



-1 and that our J(u) 



J'(u) 
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2. Finite temperature 

Using (17), the energy of a dipolar gas described by Wt is 
(after [19, 28]), ignoring factors of xj>Xz and above 

(H) T 1 1 11 
Nk B T ~ V p + 2k^ z + ^~ p + 2^ z 

(A6) 



Minimizing the energy (A6) to first order (other terms in 
the free energy do not depend on 5 or Xz/xp directly, so do 
not contribute at this level of approximation) we find Eqs. (18) 
and (19). 
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